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ABSTRACT 



We present the results of the XMM-Newton observations of five hard X-ray emitters: IGR J08262-3736, IGR J17354-3255, 
IGR J16328-4726, SAX J1818.6- 1703, and IGR J17348-2045. The first source is a confirmed supergiant high mass X-ray binary, 
the following two are candidates supergiant fast X-ray transients, SAX J1818.6— 1703 is a confirmed supergiant fast X-ray transient 
and IGR J17348-2045 is one of the still unidentified objects discovered with INTEGRAL. The XMM-Newton observations permitted 
the first detailed soft X-ray spectral and timing study of IGR J08262-3736 and provided further support in favor of the associa- 
tion of IGR J17354-3255 and IGR J16328-4726 with the supergiant fast X-ray transients. SAX J1818. 6-1703 was not detected 
by XMM-Newton, thus supporting the idea that this source reaches its lowest X-ray luminosity (=d0 32 erg/s) around apastron. For 
IGR J17348-2045 we identified for the first time the soft X-ray counterpart and proposed the association with a close-by radio object, 
suggestive of an extragalactic origin. 
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1. Introduction 



X-rays: individuals: SAX J1818. 6-1703; IGR J17348-2045; IGR J17354-3255; 



High-mass X-ray binaries (HMXBs) comprise a compact ob- 
ject orbiting a massive OB star. The compact object, usually a 
neutron star (NS), emits a conspicuous amount of X-ray radia- 
tions (up to luminosities of ~10 37 erg/s) due to the accretion of 
matter from the OB companion. Depending on the nature of the 
companion star, HMXBs can be divided into Be X-ray binaries 
(BeXBs) and supergiant X-ray binaries (SGXBs). 

In the former, the NS is in a wide eccentric orbit around a Be 
star. The X-ray luminosity is generally low (~10 32 -10 33 erg/s) 
when the NS is far away from periastron and accretes matter 
from low density regions. Remarkable increases in the X-ray lu- 
minosity (A£ x ~ 100- 1000) are displayed: (i) during the so-called 
"Type-II" X-ray outbursts, which last several orbital cycles and 
present few (if any) X-ray flux variations associated to the or- 
bital phase; (ii) at the periastron ( "Type I" outbursts), where the 
compact object is closer to the companion and accretion takes 
place through the low-velocity an d high-density wind of the Be 
star (IStella et al Jl 1 9861: iRei dl20TTh . 

In the SGXBs, the compact object moves around an early- 
type supergiant in a nearly circular orbit, and in some cases it 
is well embedded in its dense highly supersonic wind (the so- 
called "highly obscured" SGXBs). The X-ray luminosity is pow- 
ered by the accretion of the strong stellar wind onto the compact 
star and displays, on-average, less pronounced changes along the 
orbit with respect to the BeXBs. Hydrodynamic instabilities in 
the wind of the supergiant star can give rise to significant density 
gradients in the environment around the NS, leading to aperiodic 
variations of the X-ray luminosity (A^ Y ~ 10 -50) on time-scale s 
ranging from few to thousands of seconds (Negueruelal l2010h . 



The accretion of particularly dense "clumps" in the wind can 
also i nduce bright short flares that can last for a few hours (see, 
e.g jKrevkenbohm et alJ l2008. and references therein). 

A few peculiar SGXBs, sharing with the latter a remarkable 
number of similarities (orbital period, energ y spectra, propertie s 
of the companion stars; see e.g. discussion in Bo zzo et alj2 010). 
were discovered i n the late 90s (see e.g.. lYamauchi et all 1995; 
Smit h et aTll 19981) . and col lectively termed later "supergiant fast 
X-ray transients" (SFXTs; ISguera et~aLll2006h iNegueruela et al.l 
2006). At odds with the "classi cal" SGXBs, the SF XTs spend 
a large fraction of their time (iRomano et alJ l201ll) in a qui- 
escent state with a typical luminosity of 10 32 -10 33 erg/s, and 
only sporadically undergo bright outbursts (A/, x ~10 4 -10 5 ) last- 
ing a few hours and reaching peak luminosities comparable 
with those of the persistent SGXBs (Lx~10 37 erg/s; see e.g., 
Walt er & Zurita He ras 2007). The outbursts of SFXTs are asso- 
ciated to the accretion of particularly dense clumps as in other 
SGXBs, but the origin of the lower persistent luminosity and 
much more pronounced variabilit y of these sources is still a 
matter of debate din 't Zandll2005b IWalter & Zurita Herasll2007b 
Grebenev & Sunvaev 20071; Ibozzo et alj|2008d,l201 lb. 



In this paper, we report on the XMM-Newton observations 
of four HMXBs discovered with INTEGRAL and BeppoSAX. 
Among these sources, IGR J08262-3736 is classified as a classi- 
cal SGXB, IGR J 17354-3255 and IGR J 16328-4726 are candi- 
date SFXTs, and SAX J1818. 6-1703 is a confirmed SFXT. We 
also report on the XMM-Newton observation of the still unclas- 
sified INTEGRAL source IGR J17348-2045. A summary of the 
previous results for all the above sources is presented in Sect. [2] 
while our analysis and results of the XMM-Newton observations 
are presented in Sect. [3] Discussions and conclusions on the 
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properties of the soft X-ray emission from the five objects are 
reported in Sect. [4] 

2. The selected sources 

2.1. IGR J08262-3736 

IGR J08262-3736 was discovered by the hard X-ray imager 
IBIS/ISGRI on-board the INTEGRAL satel lite and reporte d for 
the first time in the fourth ISGRI catalog dBird et alj|2010t) . The 
measured fluxes of the source were 0.4 + 0.1 mCrab in the 20- 
40 keV energ y band and 0.7 ± 0.2 mCrab in the 40-100 keV 
en ergy rangdTdBird et alj|2010l) . 

iMasetti et alJ (l2010» tentatively associated the source to the 
OB-V star SS 188 located at a distance of 6.1 kpc, thus suggest- 
ing that IGR J08262-3736 was a new member of the HMXB 
class discover ed with INTEGRAL . The association was later 
confirmed by Mali zia et al.l d201 ll) using Swift/XRT data. The 
best determined XRT position is at R.A.(J2000) = 08 h 26 m 13.87 s 
and Dec.(J2000) = -37°37'1 1.03", with an associated uncer- 
tainty of 4 arcsec. The XRT spectrum could be well fit by 
using an absorbed power-law model with an absorption col- 
umn density of iVH=1.5xl0 22 cirT 2 and a power-law photon in- 
dex of T=2. The estimated flux in the 2-10 keV energy band 
was ~10" n erg/cm 2 /s, corresponding to an X-ray luminosity of 
~2.4xl0 33 erg/s. 



3 cr upper limit on its X-ray flux a factor ^30 lower than the on e 
measured in the other XRT observation (IVercellone et al.l l2009). 

The presence of these two sources was confirmed also 
thr ough the Chandra ob servation performed on 2009 February 
6 (iTomsicketalJ l2009al) . SI (=CXOU J173527.5-325554) 
was detected with an unabsorbed 0.3-10 keV flux of 
1.3xl0 _n erg/cm 2 /s, and displayed an X-ray spectrum char- 
acterized by N H ~7.5xl0 22 cm" 2 and T^O.54. The Chandra 
refined position of this source (R.A.(J2000)=17 h 35 m 27.59 s 
Dec.(J2000)=-32°55'54.4", associated 90% c.l. uncertainty 
0.64") led to the identification of the infrared counterpart 
2MASS J17352760-3255544. No catalogued optical counter- 
part was found at the Chandra position, consistent with the 
high measured extinction in the direction of the source. In 
the Chandra observation, S2 (=CXOU J173518.7-325428) 
was detected with an X-ray flux of 1.4xl0" 12 erg/cm 2 /s 
and displayed an absorbed power-law shaped spectrum with 
A^h~2.6x10 22 cm" 2 and T^O.79. The best estimated source 
position was R.A.(J2Q0Q) = 17 h 35 m 18.73 s , Dec.(72000) = 
-32°54'28.7", with an associated 90% c.l. uncertainty of 0.64". 

The more pronounced variability of SI led to the con- 
clusion that this sour ce was the most like l y counterpart 
to IGR J17354-3255 (IVercellone et al.l 120091: iTomsick et all 
20093). 



2.2. IGR J1 7354-3255 



2.3. IGR J1 6328-4726 



IGR J17354-3255 was discovered with INTEGRAL in 2006 
during a monitoring observation of the Galactic bulge region 
( Kuulkers et al. 2006), and is reported in the fourth IBIS cata- 
logue dBird et alJl2009l) and in the 54 m onths Swift/BAT hard 
X-ray catalogue dCusumano et al. I l20loh . The average X-ray 
fluxes in the 20^-0 keV and 15-150 keV energy band are 
l.lxlO" 11 erg/cm 2 /s and 2.1xl0" n erg/cm 2 /s, respectively. The 
long term monitoring of IGR J 17354-3255 carried out with the 
INTEGRAL /ISGRI and Swift/BAT suggested that the source is 
a weak persistent emitter in the hard X-rays (average flux of 
1.1 mCrab in the 18-60 keV energy band), only sporadically 
displaying relatively short flares with duration from few hours 
to ~1 day. This behaviour, together with the periodic modula- 
tion detected in the hard X-ray data at 8.4474+0.0017 days, led 
to the c onclusion that IGR J 17354-3255 is a HMXB, possibly a 
SFXT dD' Ai et al.ll201 it ISguera et al.ll201 ll) . The source is also 
positionally coincident with the high energy AGILE transient 
AGL J 1734-33 10, even though the locali zation uncertainties ar e 
still too large to claim a firm association (Bulgar elli et al.l l2009). 

In the soft X-ray domain, the source was observed twice with 
Swift /XRT a nd once with Chandra . In the most recent XRT 
observation dVercellone et al.l 120091) . two sources were found 
within the INTEGRAL error circle of IGR J17354-3255. The 
first object (SI) displayed a higher flux (l.lxlO" 11 erg/cm 2 /s, 
0.3-10 keV) and a spectrum characterized by a blackbody tem- 
perature of kT~1.4 keV and an absorption column density of 
A^h~5x10 22 cm" 2 . The second source (S2) was a factor of ~20 
fainter and no spectral information was available. In a previously 
performed XRT observation (2008 March 11), only S2 was de- 
tected at a nearly constant flux, while S 1 was not detected with a 



1 Conversions factors are lmCrab = 7.57 x 10" 12 erg/cm 2 /s in the 
2CMK) keV energy band and lmCrab = 9.42 x 10" 12 erg/cm 2 /s in the 
40-100 keV energy band, respectively. The fluxes of the source thus 
correspond to 3.0(±0.8) x 10" 12 erg/cm 2 /s and 7(±2) x 10" 12 erg/cm 2 /s, 
respectively. 



Th e source IGR J 16 328 -4726 was discovered with INTEGRAL 
by iBird et all (12007b and reported as a hard X-ray transient in 
the latest ISGRI and Swift/BAT catalogues dBird et al.l [20Tok 
ICusumano et al.l l2010h . A relatively bright outburst from the 
source was detec ted with the Swift/BAT on 2009 June 9 
(Grupe et al. 2009). The follow-up observations with the narrow 
field instrument on-board Swift, XRT, showed that the spectrum 
of the source could be well described by an absorbed power-law 
model with AT H -8xl0 22 cm" 2 and T^O.6. The estimated 0.3- 
10 keV X-ray flux was 2.4xl0" 10 erg/cm 2 /s. On that occasion, 
XRT followed the evolution o f the X-ray flux up to 4 days af- 
ter the onset of the outburst dFiocchi et al.l l2010l) . During the 
last observation performed on 2009 June 14 the source was no 
longer detected, resulting in a 3<x upp er limit on the X-r ay flux 
of 1.4xl0" 12 erg/cm 2 /s (0.3-10 keV: lFiocchi et al1l2010h . 



From the analysis of the BAT data, ICorbet et all d2010l) re- 
vealed a highly significant modulation at ~10 d. This was inter- 
preted as the orbital period of a HMXB, possibly a supergiant 
system. 

An accurate analysis of all the archival INTEGRAL observa- 
tions perf ormed in the directi on of IGR J 16328 -4726 was re- 
ported by Fioc chi et al.l d2010t) . In the hard X-ray domain (20- 
50 keV), two outbursts lasting a few hours were detected by the 
hard X-ray imager ISGRI on-board INTEGRAL. No simultane- 
ous data were available with the two X-ray monitors JEM-X. 
On both occasions the ISGRI spectrum could be well fit by us- 
ing a simple power-law model with F^2-2.6 and a flux of 2- 
3.3xl0" 10 erg/cm 2 /s (20-50 keV). The 3cr upper limit on the 
source emission outside the outbursts derived from the ISGRI 
data was 2.5xl0" 12 erg/cm 2 /s (20-50 keV). 

The fast flaring behaviour of the source, and its spectral 
properties, suggested that IGR J 16328 -4726 is a member of the 
SFXT class discovered with INTEGRAL. 
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2.4. SAX J1 81 8.6-1703 



IGRJ08262-3736 (Epic-pn) 



SAX J18 18.6-1703 is one of the confirmed SFXT sources, and 
was discovered by BeppoSAX on 1998 March 11 during a pe- 
riod of intense X-ray activity that lasted for about 2 h and 
reached a peak flux of -400 mCrab (9-25 keV, lin 'tZand etaP 
fl998h . Since then, several outbursts from the s ource were 
recor ded with INTEGRAL and Swift (see, e.g. ISidoli et aTl 
2009, and references therein). The best X-ray source position 
was provided by Chandra at R.A.(J2000)=18 h 18 m 37 s .89 and 
Dec.(J2000)=-17 °02'47.9'/(associated u ncertainty 0.6" at 90% 
confidence level; nn 'tZand et all 12 006). The estimated source 
distance is 2.1+0.1 kpc dTorreion et alJIioiOh . An in-d epth study 
of the X-ray activity of the source was car ried out by Bir d et alj 
d2009l) and IZurita Heras & Chatvl d2009l) . The first investiga- 
tion determined the best orbital period of SAX J 1 8 1 8. 6-1703 
at 30+0.1 d (in the following we also c onsider the refer - 
ence phase at 5367 1 MJD , as reported by iBird et aT] 12009). 
IZurita Heras & Chatvl J2009) also found that most of the dis- 
covered outbursts took place when the NS passed close to the 
periastron, and that the source usually remains relatively bright 
in X-rays (a few times 10~ n erg/cm 2 /s) for about ~6 d around 
this orbital phase. However, outbursts in several periastron pas- 
sages were missing. This behaviour could be reasonably well 
explained by assuming that the NS in SAX J1818.6— 1703 has 
an eccentric orbit (e=0. 3-0.4) and is sporadically accreting mass 
from the clumpy wind of its supergiant companion (see also 
Sect. [D- On 2006 October 6 SAX J1818.6-1703 was also ob- 
served with XMM-Newt on for the first time close to the apas- 
tron passage (phase 0.5 h lBird et~aDl2009h . but was not detected. 
iBozzo eta! (2008a) determined a 3<x upper limit on its unab- 
sorbed X-ray flux of l.lxl0~ 13 erg /cm 2 / s (the effective expo- 
sure time was ~13 ks). 



2.5. IGR J1 7348-2045 

IGR J17348-2045 is an unclassified sourc e discovered wit h 
INTEGRAL and reported for the first time bv lBird etall d2007l) . 
The best determined position of the source with ISGRI is at 
R.A.(J2QQQ) = 17 h 34 m 57 s , Dec.(J2QQQ) = -20j^4J^with a 
90% c.l. associated uncertainty of 3.6 arcmin dBird et al. 2010). 
The estimated fluxes in the 20^40 keV and 40-100 keV energy 
bands are 0.3 + 0.1 mCrab and 0.9 ±0.1 mCrab, respectively. 



3. Data analysis and results 

For the present study we used all the available XMM-Newton 
observations of the five sources obtained through the guest ob- 
servation programs awarded to our research group. 

XMM-Newton observation data files (ODFs) were processed 
to produce calibrated event lists using the standard XMM- 
Newton Science Analysis System (v. 11.0). We used the epproc 
and emproc tasks to produce cleaned event files from the EPIC- 
pn and MOS cameras, respectively. EPIC-pn and EPIC-MOS 
event files were filtered in the 0.5-12 keV and 0.5-10 keV en- 
ergy range, respectively, to exclude high background time in- 
tervals. The effective exposure time for each observation and 
camera is specified for each source in the following sections. 
Lightcurves and spectra for the five sources and the relative 
backgrounds were extracted by using regions in the same CCD. 
The difference in extraction areas between source and back- 
ground was accounted for by using the SAS backscale task for 
the spectra and lccorr for the lightcurves. All EPIC spectra were 
rebinned before fitting in order to have at least 25 counts per bin 




5000 



10 4 1.5X10 4 2xl0 4 2.5x10" 

Time (s) 



Fig. 1. XMM-Newton background-subtracted lightcurve of 
IGR J08262-3736 in two energy bands and the corresponding 
hardness ratio. The latter is calculated as the ratio between the 
source count rate in the 2-12 keV and 0.5-2 keV energy bands. 
The time bin is 150 s. 



and, at the same time, to prevent oversampling of the energy 
resolution by more than a factor of three. Where required, we 
barycenter-corrected the photon arrival times in the EPIC event 
files with the barycen tool. Throughout this paper, uncertainties 
are given at 90% c.l., unless stated otherwise. 

3.1. IGR J08262-3736 

XMM-Newton observed the field-of-view (FOV) around 
IGR J08262-3736 on 2010 October 16 for a total exposure time 
of ~25 ks. All three EPIC cameras were operated in full frame 
mode. A bright point source was detected at a position (obtained 
from the edetect.chain tool) RA.(J2000)=08 h 26 m 13 s :7 and 
Dec.(J2000)=-37°37'll'.'58 (the nominal positional accuracy of 
the EPIC cameras is assumed her eafter to be ~2 arc sec unless 
otherwise stated; see discussion in iPavan et alJl201ll) . compati- 
ble with that of the soft X-ray cou nterpart of IGR J08262-3736 
identified by [Malizia et alj d201 ll) . On the EPIC-pn, the source 
lied close to the border between two CCDs and thus we chose 
the extraction region in order to avoid any contribution from the 
CCD rirrQ. The source lightcurve in two energy bands is shown 
in Fig.[T] together with the hardness ratio. The source displayed 
a moderate variability, with two relatively small flares occurring 
about 1.5xl0 4 s after the beginning of the observation. The 
hardness ratio and the hardness intensity diagram of the source 
(see Fig. Q]and[2| did not show clear evidence of variations in 
the source spectral properties with the count-rate. 

An acceptable fit to the EPIC-pn spectrum of the source 
accumulated over the entire exposure time available could 
not be obtained using a simple absorbed power-law model 
(Y 2 ed /d.o.f. = 1.7/1 34). The fit was significantly improved by 
adding a blackbody component at the lower energies (here- 
after, BB) or by the introduction of a partial covering model. 
In order to constrain better the spectral parameters, we also ex- 
tracted the spectra from the EPIC-MOS cameras and the aver- 
age INTEGRAL /ISGRI spectrum of IGR J08262-3736 (total ef- 



2 At the time of the XMM -Newton observatio n, the refined XRT po- 
sition of IGR J08262-3736 dMalizia et al.ll201ll) was not available. 
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Table 1. Results of the spectral fit parameters obtained from the average EPIC+ISGRI spectra of IGR J08262-3736 (see Sect. [XT). 



Model Na 


r 




/ 


kT BB 


^BB 


F 0.5-2 keV 


^2-10 keV 


^20-40 keV 


Cmosi (CmOS2) 


ClSGRl 




1 1.4±0.1 
0.60±0.07 

2 0.80±0.06 


1.52±0.05 

1.3±0.1 

1.8±0.1 


3.8±0.7 


0.63±0.04 


0.10±0.01 
1.2±0.2 


132+ 11 

J -33 

0.18±0.03 


16.0 
1.0 

4.5 


6.0 

5.3 
6.5 


4.0 

3.6 
4.4 


1.03±0.03 (1.01±0.03) 
1.00±0.03 (1.00±0.03) 
1.00±0.03 (1.00±0.03) 


0.7±0.2 

w& 

1.2±0.4 


1.1/380 
0.9/380 
0.9/380 



Notes. The leftmost column indicates the spectral model used in Xspec. Model 1 is phabs*(BB+pow), model 2 is phabs*pcfabs*pow (phabs is the 
absorption component, pow the powerlaw, BB the blackbody, and pcfabs the partial covering). We indicated with T the power-law photon index, 
/Vh the absorption column density in units of 10 22 citT 2 , / the covering fraction in the partial covering model (the local absorption component is 
indicated with N Hb in units of 10 22 cm -2 ), and kT BB (R B b) the temperature (radius) of the BB component in keV (km). For R BB we considered a 
source distance of 6.1 kpc (see Sect.fT}. For the model that includes the BB component (Model 1), we report both the best fit values obtained with 
a hotter (kT~l keV) and colder (kT~0. 1 keV) emitting surface. We also reported the normalization constants Cmosi, Cmos2> and Qsgri for the 
MOS1, MOS2 and ISGRI spectra, respectively (we fixed the EPIC-pn constant to 1 as a reference). The unabsorbed fluxes (units of 10~ 12 erg/cm 2 /s) 
in the 0.5-2 keV, 2-10 keV, and 20-40 keV are indicated in all cases with F 5 _ 2 keV, ^2-10 keV, and ^20-40 keV, respectively. 



IGR J08262-3736 (Epic-pn) 




0.5 1 2 

Intensity [cts/s] 



Fig. 2. Hardness-intensity diagram of IGR J08262-3736 ob- 
tained by using EPIC-pn data. For this plot, data were rebinned 
in order to have S/N>5 in each bin. 



IGRJ08262-3736 (EPIC and ISGRI) 




1 1 

Energy (keV) 



Fig. 3. Average XMM-Newton /EPIC and INTEGRAL /ISGRI 
spectra of IGR J08262-3736. The best-fit model (solid line) is 
obtained with a partial-covering power-law model. The residuals 
from the fit are shown in the bottom panel. 



fective exposure time 1.2 Ms) from the Heavens on-line toofl. 
The combined EPIC+ISGRI spectra of the source are shown in 
Fig. [3] The results of all the fits are reported in Table Q] A nor- 
malization constant with respect to the EPIC-pn was included in 
all the fits to take into account the inter-calibration between the 
instruments and the possible variability of the source. No signif- 
icant evidence for coherent modulations was found in the EPIC 
data of IGR J08262-3736. 



3.2. IGR J1 7354-3255 

The FOV around IGR J17354-3255 was observed with XMM- 
Newton on 201 1 March 6 for a total exposure time of 20 ks. The 
EPIC-MOS 1 and EPIC-pn cameras were operated in full-frame, 
while the EPIC-MOS2 was in small-window mode. No signifi- 
cant time intervals were affected by a high flaring background, 
and thus in the following analysis we retained the entire expo- 
sure time available for the three instruments. 

Among the two candidate soft X-ray counterparts of 
IGR J 17354-3255 identified previously with Chandra (see 
Sect. 12.21 1. SI was not detected by the three EPIC cameras. We 
calculated an upper limit on the emission from the source by 
using the ximage tool uplimit, and a circular extraction region 
of 15 arcsec centered on the Chandra position of the source (see 
Sect. 12.21 ). We obtained a 3 cr upper limit on the source count-rate 
of 0.002 cts/sec in the 0.5-10 keV energy band (effective expo- 
sure time 19.0 ks). This corresponds to an observed X-ray flux 
of 7xl0~ 14 erg/cm 2 /s(we assumed a spectral index of F=1.7 and 
an absorption colum n density of A^h=7x10 22 cirT 2 as reported 
bv lD'AietalhollI) . 

S2 was outside the EPIC-MOS2 FOV, but clearly de- 
tected by the EPIC-pn and the MOS1. In the former instru- 
ment the source was located on the gaps between two CCDs, 
and thus only data from the EPIC-MOS 1 are considered in 
the following analysis. The best determined source position 
with the edetect.chajn tool is R.A.(J2000)=17 h 35 m 18 s .48 and 
Dec.(J2000)=-32°54'30'.'2, compatible with that determined pre- 
viously with Chandra (see Sect. I2.2| FI. We extracted the spec- 
trum of S2 from a cir cular region centered on the Chandra 
position of the source (Qbrrisick et al] |2009al) . Due to the lim- 
ited statistics of the data, the spectrum was grouped to have at 



3 http://www.isdc.unige.ch/heavens. 

4 The source is relatively faint for the MOS1, and a position un- 
certainty as large as 4"can be expected in this case (see discussion in 
Pavanetal 1l201lh . 
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least 5 counts per brrQ and fit with the C-statistics (ICashll 19791) 
by using an absorbed power-law model. We fixed the absorp- 
tion column density to the Galactic value expected in the di- 
rection of the source (jV ff ^1.4xl0 22 cm~ 2 . lKalberla et al.ll2005t 
iDickev & Lockmanl 1990) and measured a power-law photon in- 
dex of r=0.9+0.3 (C-statistics/d.o.f.=29.4/36). The correspond- 
ing 0.5-10 keV observed flux was 2.4xl0~ 13 erg/cm 2 /s. No sig- 
nificant evidence for coherent modulations was found in the 
EPIC data. 

3.3. IGR J1 6328-4726 

IGR J16328-4726 was observed by XMM-Newton on 2011 
February 20 for a total exposure time of ~22 ks. Filtering 
for the flaring background intervals resulted in an effective 
exposure time of 14.7 ks for the EPIC-pn and 21 ks for 
the two MOS cameras. The best determined source position 
with the edetect_chain tool is R.A.(J2000)=16 h 32 m 37.68 and 
Dec.(J2000)=-47°23'3 9.48", compatible with that determined 
previously with Swift dGrupe et al.l l2009h . The lightcurves ex- 
tracted from the three instruments and corrected for the back- 
ground are shown in Fig. [4] The source displayed a clear vari- 
ability, with relatively small flares (A^ x <10) occurring during 
periods of low-level X-ray activity. The average EPIC-pn spec- 
trum extracted during the observation is shown in Fig. [5] The 
best fit to this spectrum (y 2 ed /d.o.f. = l. 06/105) is obtained by us- 
ing an absorbed power-law model. The measured absorption col- 
umn density is AfH=(17.5±1.0)xl0 22 crrT 2 , and the power-law 
photon index is r=1.5±0.1. The average unabsorbed 2-10 keV 
flux is 1.7xl0~ n erg/cm 2 /s. We searched for possible spectral 
variations during the rise and decay of the small flares by using 
the hardness-intensity diagram (HID) of the source (see Fig. [6]). 
The Spearman rank correlation coefficient calculated for the 97 
data points comprised in the plot is 0.39, thus indicating a statis- 
tically significant correlation between the HR and the source in- 
tensity (the corresponding null hypothesis probability is <10~ 4 ). 
We carried out a count-rate resolved spectral analysis to check 
this result, but did not find any significant variation of the spec- 
tral parameters (A^h, F) with the source intensity. This can be as- 
cribed also to the relatively low statistics of the data. We checked 
that this result could not be improved by fitting simultaneously 
pn and MOS data extracted during the same time intervals (we 
used the good time intervals, GTIs, of the EPIC-pn to extract si- 
multaneous spectra from the two MOS cameras). The results of 
this analysis are shown in Fig. [7] We report there only the results 
corresponding to one particular selection of the source high and 
low count-rate time intervals. Different choices of these intervals 
did not lead to significant changes of the results. 

As shown in Fig. [4] the two EPIC-MOS cameras began ob- 
serving the source about 4 ks earlier than the EPIC-pn. Beside 
the small flares visible from all the instruments, the EPIC-MOS 
showed an intriguing drop of the source count-rate ~3 ks after 
the beginning of the observation. We performed a spectral analy- 
sis of the combined MOS 1 and MOS2 data extracted during this 
evenfl but the relatively low statistics did not permit to carry out 
any detailed spectral analysis (only 60 counts were recorded dur- 
ing the -0.7 ks corresponding to the event). The derived spectral 
parameters were fully consistent with those measured from the 



5 We performed this minimal grouping to avoid issues related to 
empty channels, see heasarc.nasa.gov/xanadu/xspec/XspecManual.pdf. 
A check was carried out a posteriori to verify that this did not affect the 
results. 

6 See http://xmm.esac.esa.int/sas/current/documentation/threads/ 
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Fig. 4. XMM-Newton background-subtracted lightcurve of 
IGR J16328-4726. Data are shown in the 0.5-12 keV and 0.5- 
10 keV energy range for the EPIC-pn and EPIC-MOS, respec- 
tively. The time bin in all cases is 100 s. 
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Fig. 5. Average XMM-Newton spectrum of IGR J 16328 -4726 
extracted from the EPIC-pn data. The best-fit model is obtained 
with an absorbed power law (see text for details). The residuals 
from the fit are shown in the bottom panel. 



average spectrum, in turns compatible (to within the errors) with 
those already obtained from the spectral fit to the EPIC-pn data. 
The highest and lowest unabsorbed fluxes of IGR J 16328 -4726 
recorded by the EPIC cameras were 5.2xl0~ n erg/cm 2 /s and 
6.4xl0~ 12 erg/cm 2 /s, respectively. No significant evidence for 
coherent modulations was found in the EPIC data of the source. 



3.4. SAX J1 81 8.6- 1703 

The XMM-Newton observation of SAX J1818.6-1703 reported 
here was performed on 2010 March 21 starting from 12:20 UTC 
for a total exposure time of 45 ks. The EPIC-pn camera was 
operated in Full Frame mode, while the EPIC-MOS 1 and EPIC- 
MOS2 were operated in Small Window and Fast Uncompressed 
mode, respectively. The source was not detected by the EPIC 
cameras. To estimate an upper limit on the source X-ray flux, we 
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Fig. 6. Hardness-intensity diagram of IGR J16328-4726 real- 
ized by using EPIC-pn data. For this plot, data were rebinned in 
order to obtain in each bin a S/N>5. 




Fig. 8. XMM-Newton image of the FOV around 
IGR J17348-2045 (EPIC-MOS image in the 0.5- 
10 keV energy band). The two soft X-ray counterparts of 
IGR J17348-2045 detected within the INTEGRAL error cir- 
cle of the source are indicated with blue circles. The radio 
object NVSSJ173459-204533, spatially coincident with 
XMMU J173458. 8-204530, is also indicated with a blue cross. 



IGRJ 16328-4726 (EPIC spectra) 




Energy (keV) 

Fig. 7. EPIC-pn and EPIC-MOS spectra of IGR J 16328 -4726. 
The higher three spectra were extracted by selecting only the 
time intervals in which the EPIC-pn source count-rate was 
>1.5 cts/s (0.5-12keV, see Fig. @J. The lower three spectra 
were extracted during time intervals in which the count-rate 
was <0.7 cts/s (same energy band). GTIs from the EPIC-pn 
were used also to extract the corresponding EPIC-MOS spec- 
tra. The best-fit model (solid line; ;^ 2 ed /d.o.f.=0. 9/332) is ob- 
tained with an absorbed power-law (the residuals from the fit 
are shown in the bottom panel). We obtained for the high- 
est spectra N H =(17.4±1.7)xl0 22 cirT 2 and r=1.6+0.2; for 
the lower spectra we obtained Nh=(18.2+1.8)x10 22 crrT 2 and 
T= 1.5+0. 2. The corresponding observed 2-10 keV X-ray fluxes 
were 2.2xl0~ n erg/cm 2 /s and 4xl0~ 12 erg/cm 2 /s, respectively. 



first filtered the EPIC data for the high flaring background and 
then applied the same technique as described in Sect. 13.21 As 
the observation suffered from a very high level of flaring back- 
ground, the effective exposure time available for the EPIC-pn 
was 4.0 ks. We estimated from these data a 3 cr upper limit on the 
source X-ray count-rate of 0.0052 cts/s, and converted this value 
into a 0.5-10 keV unabsorbed flux of 3.0xl0~ 13 erg/cm 2 /s by 



assuming the same spectral model used in lBozzo et al.l d2008b). 
This allows an easier comparison with the previous XMM- 
Newton non-detection of the source (see Sect. 12.41 ). We also 
checked that a less conservative selection of GTIs for the EPIC- 
pn data would not affect significantly the derived upper limit on 
the X-ray flux of SAX J18 18.6- 1703. 

According to the ephemeris given by iBird et aD (|2009), 
the epochs of the present XMM-Newton observation of 
SAXJ1818.6-1703 corresponds to phase 0.53±0.18. Even 
though the associated uncertainty is rather large, this result sug- 
gests that the XMM-Newton observation took place around the 
system apastron. The upper limit on the flux is indeed compa- 
rable with that estim ated by the previou s non-detection of the 
source at phase 0.5 1 dBozzo et 



le previous 
aDl2008b). 



3.5. IGR J1 7348-2045 

IGR J 17348 -2045 was observed by XMM-Newton on 2011 
March 2 for a total exposure time of 20 ksec. The three EPIC 
detectors were operated in Full Frame with medium filter. 

The observation was affected by a high flaring background, 
resulting in an effective exposure time of 10 ks for the two 
MOS cameras and 4 ks for pn. We searched simultaneously 
for the detected sources in all the cleaned EPIC images by 
using the edetect_chain SAS tool. Inside the INTEGRAL 
error circle around IGR J 17348 -2045 we found two soft X-ray 
sources located at RA = 17 h 34 m 58.80 s , Dec = -20°45'30.96", 
and RA = 17 h 34 m 49.20 s , Dec = -20°42'44.64", respec- 
tively. According to the XMM-Newton convention, we 
named the two sources XMMU J173458. 8-204530 and 
XMMUJ173449.2-204244. The former source is clearly 
visible in all the images of the individual instruments and 
located close to the center of the INTEGRAL error-circle; 
the latter is detected only when all EPIC images are used 
simultaneously and located at the very rim of the INTEGRAL 
error circle (see Fig. [8]). For this latter source no meaningful 
spectral information could be extracted, and we thus assume 
in the following that XMMU J173458.8-204530 is the true 
counterpart to IGR J 17348 -2045. For this source, the statistic 
was enough to perform a spectral analysis. We extracted the 
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IGRJ 17348-2045 (Epic-MOS and ISGRI) 
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Fig. 9. EPIC-MOS and ISGRI spectra of 
XMMUJ173458. 8-204530, the most likely counterpart to 
IGR J 17348 -2045. The best fit is obtained by using an absorbed 
power-law model (see Sect. l3.5l for details). 



K=9.019 mag), respectively. These two objects are located at 
the rim of the XMM-Newton error circle. In addition, the large 
absorption measured in the X-ray domain, whi ch corresponds 
to an optical extinction of ^8 mag according to iGiiver & Ozell 
(2009), suggests that an optical counterpart might not be 
expected. Further observations are needed to confirm these 
associations (e.g. with Chandra ). Fro m the 1 .4 GHz conti nuum 
NRAO Very Large Sky Survey (NVSS ICondon et al.|[l998l) . we 
also noticed that the object NVSS J173459-204533 lies rela- 
tively close to the XMM source (~4.8"). The best determined 
position of the radio object is at RA = 17 h 34 m 59.08 s (±0.09 s), 
Dec = -20°45'33.8" (±1.2 arcsec), and the integrate d flux 
density at 1.4 GHz is 13.2 ± 0.6 mJy (ICondon et alJll998h . 



4. Discussion and conclusions 

We reported here on the observations of five hard X-ray 
emitters obtained with XMM-Newton. Among these, four 
(IGR J08262-3736, IGR J17354-3255, IGR J16328-4726, 
SAXJ1818.6-1703) are HMXBs, whereas IGR J17348-2045 
is still unidentified. 



three MOS and pn spectra using a circular extraction region 
centered on the best determined source position (see above). 
The background spectrum was extracted from the closest 
source-free region to XMMU J173458. 8-204530 (we checked 
that different choices of the background extraction region did 
not affect significantly the spectral results). On the EPIC-pn 
the source was located on the gap between two CCDs, and 
thus no spectrum could be extracted. A simultaneous fit to 
the EPIC-MOS spectra with an absorbed power-law model 
gave /Vh=17^x10 22 cirr 2 , F=1.5+^ <^ 2 ed /d.o.f.= 0.7/16). 
The estimated absorption column density is much higher than 
the Galactic value expected in the direct ion of the source 
(~2xl0 21 cm" 2 ; iDickev & Lockmanl [l990h . The unabsorbed 
2-10 keV flux is F2-10 kev=2.2xl0 -12 erg/cm 2 /s. No significant 
evidence for coherent modulations was found in the EPIC data 
of IGR J17348-2045. For this source we also retrieved the 
long-term ISGRI spectrum from the Heavens on-line tool (effec- 
tive exposure time 3.9 Ms) and performed a combined fit with 
the MOS spectra (see Fig. [9]). A normalization constant with 
respect to MOS was included in the fit to account for the inter- 
calibration between the instruments and possible variability of 
the source. The best fit model (y 2 ed /d.o.f.= 0.7/18) was obtained 
using the same absorbed power-law model as described above, 
with values of the absorption column density and power-law 
photon index fully compatible to within the uncertainties (no 
significant improvement on the uncertainties of the spectral 
parameters could be obtained). The normalization constant is 
0.4 ^9-?, thus supporting the idea that the source is a persistent 
hard X-ray emitter (we verified that fixing the normalization 
constant to 1 for both the ISGRI and MOS spectra would not 
affect significantly the results of the fit). The estimated 20- 
40 keV and 40-100 keV X-ray fluxes were 2xl0~ 12 erg/cm 2 /s 
and 5.2xl0~ 12 erg/cm 2 /s, respectively. The se are compatible 
with the values reported bv lBird et alJ d2009l) . We also searched 
for counterparts to XMMU J173458.8-204530 in the optical, 
infrared and radio domain. The closest catalogued objects 
to the XMM source in the optical and infrared domains are 
USNO-B 1.0 0692-0466331 (Bl = 15.82 mag, Rl = 12.97 mag, 
B2=15.40 mag, R2=14.35 mag, 1=11.70 mag) and 
2MASSJ17345863-2045292 (J=10.532 mag, H=9.446 mag, 



4.1. IGR J08262-3736 

The XMM-Newton observation of this source revealed a timing 
and spectral behaviour that are typical of classical SGXBs. The 
EPIC-pn lightcurve shows a relatively moderate variability on 
time scales of hundreds seconds with variations in the X-ray 
flux of a factor ~3-4. These variations are usually ascribed to 
changes in the mass accretion rate due to instabilities and den- 
sity fluctuations in the supergiant wind (see Sect. [TJ. The EPIC 
spectra could not be fit with a simple absorbed power-law model, 
and clearly showed the presence of an excess in the residuals 
below 2-3 keV. The presence of this "soft excess", requiring 
an additional spectral component to obtain an acceptable fit to 
the data, is believed to be a ubiquitous feature of binary sys- 
tems hosting accreting NSs dHickox et alj |2004). Its detectabil- 
ity is mostly related to the absorption column density in the 
direction of the source and its flux. In wind accreting binaries 
with X-ray luminosities comparable to that of IGR J08262-3736 
(^3xl0 34 erg/s, assuming a distance of 6.1 kpc), the soft compo- 
nent is most likely originated by thermal X-ray photons close to 
the NS surface or by pho toionized or collisio nally heated diffuse 
gas in the binary system (Hick ox et alJ l2004). 

The statistics of the EPIC spectra below 2-3 keV shown in 
Sect. 13.11 were unfortunately too low to clearly distinguish be- 
tween these possibilities. An acceptable fit could be obtained by 
using either a hot (kT^l.2 keV) compact (^0.1 km) BB compo- 
nent, or a colder {kT ^0.1 keV) and more extended (~ 130 km) 
thermal emission. The latter model would probably be more 
plausible in the case of IGR J08262-3736 as low luminosity 
wind accreting systems are not expected to d isplay hot emitting 
spots on the NS surface (see discussion in iBozzo et al.l 12010. 
and references therein). As an alternative interpretation, we also 
showed in Table Q] that a partial-covering model would provide 
an acceptable description of the X-ray spectrum. In this case, the 
soft excess would be produced by the effect of partial obscura- 
tion of the emission from the NS by the surrounding high density 
material (see e.g., Tomsick et al. 2009b). 

The combined fit of the EPIC spectra with that obtained from 
the long-term monitoring of the source with ISGRI, revealed that 
the models described above can also account for the higher en- 
ergy emission recorded from the source (up to ~60 keV). The 



7 



Bozzo et al.: XMM-Newton observations of four HMXBs and IGR J17348-2045 



normalization constants between all the instruments were com- 
patible with unity, thus suggesting that IGR J08262-3736 is 
characterized by a virtually constant persistent X-ray flux. This 
is expected for a classical SGXB. 

4.2. IGR J1 7354-3255 

The XMM-Newton observation detected only one of the two 
possible counterparts to the INTEGRAL source identified be- 
fore with Swift and Chandra (see Sect. I3.21 >. For the most likely 
counterpart (SI), the non-detection by XMM-Newton sets a 3<x 
upper limit of 7xl(T 14 erg/cm 2 /s on its observed 0.5-10 keV 
X-ray flux. This is a factor of ~6 lower than the previously de- 
termined upper limit (see Sect. I3.2l i. This observation thus in- 
creases significantly the dynamic range in the X-ray luminosity 
of IGR J17354-3255, supporting the classification of this source 
as a SFXT. 

We also note that the X-ray flux of the source S2 measured 
by XMM- Newton is roughly comp atible with that reported pre- 
viously bv lVercellone et a l. (2009), thus supporting the idea that 
this is a persistent object unrelated to the INTEGRAL source. 

4.3. IGR J1 6328-4726 

The XMM-Newton observation of this source revealed a vari- 
ability on time scales of hundreds of seconds (see Fig. |4|l 
that is very reminiscent of that observed from the SFXT pro- 
totyp es IGR J 1739 1-3021 and IGR J08408 -4503 (IBozzo et al l 
2010). This provides further support to the association of 
IGR J16328-4726 with this class of objects. 

At odds with some other SFXTs observed in quies- 
cence, we could not detect any soft excess in the X-ray 
spectrum of the source. According to the discussion in 
Sect. 14.11 this spectral component might have gone undetected 
in IGR J16328-4726 due to its relatively high absorption col- 
umn density (YVh~1-8x10 23 crrT 2 ) and the lower exposure time 
with respect to the XMM-Newton observations of other quies- 
cent S FXT sources with a similar X-ray luminosity (IBozzo et al.l 
120101) . 

As a final remark we note that the lowest X-ray flux mea- 
sured by XMM-Newton is compatible with the upper lim it re- 
ported previously with Swift /XRT bv lFiocchi et alJ ([20 10), thus 
suggesting that this is the real quiescent emission level of the 
source. 

4.4. SAX J1 81 8.6-1 703 

This source is one of the confirmed SFXTs, and has been ob- 
served in outburst several times with INTEGRAL and Swift. The 
properties of its soft X-ray emission (<10 keV) during outbursts 
were investigated in detail thanks to the observations performed 
with the XRT on-board Swift, but so far no observations with 
a focusing X-ray telescope endowed with a high sensitivity was 
able to detect the source in quiescence. The XMM-Newton obser- 
vation reported here and the one performed previously in 2006 
(IBozzo et a l. 2008b) were carried out close to the system apas- 
tron and provided a similar 3<x upper limit on the source X-ray 
luminosity of ^2xl0 32 erg/s (assuming a distance of 2.5 kpc). 

The nature of the prolonged quiescent state of 
SAX J 1 8 1 8. 6-1703 around the apastron is presently un- 
known, due to the lack of proper spectral information. As 
the orbital period of the system is ~30 days, the presence of 
X-ray eclipses during apastron would imply relatively strong 



constraints on the system inclination. A high eccentricity (e) 
might help to reduce the accretion rate at the apastron, but 
the value of e inferred for SAX J1818.6-1703 (0.3-0.4) can 
hardly explain the large dynamic range in the X-ray luminosity 
displa yed by the source (Bird et al. 200*91; IZurita Heras & Chatvl 
2009). A further reduction of the mass accretion rate would 
thus require some additional mechanism to be at work at 
this orbital phase. Proposed models involve the inhibition 
of the accretion by the NS centrif u gal and/or the mag netic 
barrier dGrebenev & Sunvaevl 120071: Bozzo et al.l l2008bl) . or 
the presence of a highly structured wind from the companion 
star which is extremely rarefied around apastron and permits 
at this phase only accretio n at a very low-level (see, e.g. 
IZurita Heras & Chatvl |2009, and references therein). Further 
observations of SAX J 1 8 1 8.6- 1 703 around orbital phase 0.4-0.6 
with high sensitivity X-ray telescopes, like XMM-Newton and 
Chandra, are required to clarify the origin of the lowest X-ray 
luminosity state reached by this source. 

4.5. IGR J1 7348-2045 

This source was never observed before in the soft X-ray do- 
main. The XMM-Newton data reported here allowed us to iden- 
tify the counterpart to the INTEGRAL source and revealed 
that the source X-ray emission is intrinsically highly absorbed 
(^10 23 cirT 2 , see Sect. \3.5\ and well described by a power-law 
model with F~1.5. A simultaneous fit with the long-term ISGRI 
spectrum also supports the idea that the source is a persistent 
hard X-ray emitter. 

The uncertainty in the EPIC-MOS position does not allow 
to firmly identify the counterparts at other wavelengths, but 
we remark here the possible association with the radio object 
NVSSJ173459-204533. If this association will be confirmed 
by future observations performed with higher spatial-resolution 
instruments, we propose that IGR J17348-2045 could be a 
highly absorbed active Galactic nuclei (AGN). A number of 
thes e sources were re cently discovered with INTEGRAL (see, 
e.g. iRicci et "aIl l20Tll and references therein). Further multi- 
wavelength observations of this source are required in order to 
establish its real nature. 
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